Signals of Warped Extra Dimensions at the LHC 
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Abstract. We discuss the signatures of the spin-2 graviton excitations predicted by the Randall- 
Sundrum model with one warped extra dimension, in dilepton and diphoton production at LHC. 
By using a specific angular analysis, we assess the ranges in mass and coupling constant where 
such gravitons can be discriminated against competitor spin-1 and spin-0 objects, that potentially 
could manifest themselves in these processes with the same mass and rate of events. Depending 
on the value of the coupling constant to quarks and leptons, the numerical results indicate graviton 
identification mass ranges up to 1.1-2.4 TeV and 1.6-3.2 TeV for LHC nominal energy of 14 TeV 
and time-integrated luminosity of 10 and 100 fb , respectively. 
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INTRODUCTION 

Originally, extra spatial dimensions were proposed to address the mass scale hierarchies 
affecting the Standard Model (SM), and requiring parameter fine tuning, in particular the 
so-called gauge hierarchy M^w <^ Mpi ~ 10 16 TeV. These scenarios predict the existence 
of heavy new particles, or excitations of the SM particles, that could be revealed as 
narrow peaks in cross sections measured at the LHC if their masses were within the 
experimental kinematical reach. Observability will depend on the values of these heavy 
'resonance' masses and coupling constants to the SM matter, that assure sufficiently 
high event rates. 

Actually, in this regard, the discovery reach on an individual scenario is the maxi- 
mum value of the corresponding resonance mass, Mr, for which the peak can be ob- 
served. However, different non-standard models can give, for appropriate values of their 
parameters, peaks at the same M = Mr and the same number of events under the peak. 
Accordingly, one defines for any non-standard scenario an identification mass range (of 
course, included within the discovery range), where it can be unambiguously discrimi- 
nated as the one really underlying the peak, by excluding the other models as potential 
sources of that same peak. The determination of the spin of a discovered heavy reso- 
nance is, therefore, crucial for its identification. We will here discuss the identification 
of the spin-2 graviton excitation predicted by the Randall-Sundrum (RS) model with one 
warped extra dimension [1], against the spin-1 and the spin-0 hypotheses for a heavy 
neutral resonance discovered in dilepton and diphoton inclusive production at the LHC: 

p + p -> l + l~ +X (l = e,n) and p + p -+yy+X. (1) 



We will model the spin-1 hypothesis by the Z's predicted by extended electroweak 
gauge symmetries [2], and the spin-0 one by the sneutrinos (v) envisaged by 7?-parity 
violating SUSY extensions of the SM [3]. 'Confusion' domains of their respective 
parameter spaces allowed by current experimental constraints exist, in which '^-channel' 
exchanges of the above mentioned particles can in the process (1) produce narrow 
peaks in the dilepton invariant mass with the same values of mass and event rates. 
Therefore, for the spin-2 hypothesis discrimination against the two alternative ones, 
more detailed information must be supplied, i.e., that embodied in the characteristic 
angular distributions of the different scenarios. 

BASIC OBSERVABLES AND ANGULAR ANALYSIS 

The basic observable for a heavy resonance discovery at an invariant dilepton (or dipho- 
ton) mass M = M R (with in our case R = G,Z' ', v denoting graviton, Z' and sneutrino, 
respectively) is the production cross section governing the rate of events 

/Zcut rM R +AM/2 ry m . dx 

dz m dy , (2) 

-zcm J m r - am /2 Jy min dMdydz 

and the differential angular distribution 

d(7 i-Mr+AM/2 ry m - dx ^(7 

— =/ dM I dy. (3) 

dz J Mr- am ji Jy min dMdydz 

In Eqs. (2) and (3), z = cos cm and y define the lepton-quark (or photon-quark) angle 
in the dilepton (or diphoton) center-of-mass and the dilepton rapidity, respectively, and 
cuts on phase space due to detector acceptance are indicated. For integration over the full 
phase space, the limits would be z C ut = 1 and y max = — y m in = ^og(y/s/M) with y/s the 
LHC collider center-of-mass energy. Furthermore, AM is an invariant mass bin around 
Mr, reflecting the detector energy resolution, see for instance Ref. [4]. To evaluate the 
number Ns of resonant signal events, these equations must be multiplied by the time- 
integrated luminosity, for which we will take 100 and 10 fb _1 , and by the foreseen 
reconstruction efficiencies (90% for both electrons and muons [5]). Also, the typical 
experimental cuts on rapidities and transverse momenta must be applied (p± > 20 GeV 
and |rj | < 2.5 for both leptons). Finally, with Nb the number of 'background' events in 
the AM bin, determined by the SM predictions, the criterion Ns = 5y/N~B or 10 events, 
whichever is larger, will be adopted as the minimum signal for the peak discovery. Of 
course, the determination of discovery and identification reaches on the different non- 
standard models requires the expression of Eqs. (2) and (3) in terms of convolutions of 
the pertinent partonic cross sections times parton distribution functions, and for the latter 
we will choose the CTEQ6 ones of Ref. [6]. 

Since the z-dependence of the cross sections for graviton, Z' and sneutrino exchanges 
are different, the discrimination of spin-2 from spin-1 and spin-0 might be attempted by 
the 'direct' comparison of Eq. (3) for the three hypotheses [7, 8]. 

In practice, due to the completely symmetric pp initial state at the LHC, the deter- 
mination of the sign of z from the measured events may not be fully unambiguous. An 



observable potentially avoiding this ambiguity, and which we adopt to perform the an- 
gular analysis, is the z-evenly integrated center-edge angular asymmetry [9, 10]: 



A CE = — with Oce = \ f - ( [ * + H)} ^-dz. (4) 

° U-Z* \J-Zcat JZ* J] aZ 

In Eq. (4), < z* < z cu t defines the separation between the "center" and the "edge" 
angular regions, and is a priori arbitrary. However, its actual value can be "optimized" 
in the numerical analysis. A further potential advantage of Ace is that, as consisting of 
ratios of integrated cross sections, it could be less sensitive to systematic uncertainties, 
such as those stemming from particular choices of parton distributions, ^-factor values, 
etc.. 1 



NON-STANDARD MODELS AND ANGULAR DISTRIBUTIONS 
RS model with one compactified warped extra dimension 

This simplest version consists of one warped extra dimension, y, two three-dimensional 
branes placed at a compactification relative distance kR c in y, and the specific 5-D metric 
[1] 

ds 2 = exp(-2&|y|) ri^ v dx^dx v -dy 2 , (5) 

where r\^ v is the usual Minkowski tensor and k > is a dimensionful constant. The 
SM fields are assumed to be localized to the so-called TeV brane, while gravity can 
propagate in the full 5-D space, in particular on the other brane, the Planck brane, 
in which the effective 4-D mass scale is Af H = \/^%kG^ = 2.44 x 10 15 TeV. With 
M* the 5-D mass scale, analogously related to the cubic root of the 5-D gravitational 

constant, Einstein's equations imply the relation M P1 = (M*/k) ( 1 — exp ( —2knR c ) ) , and 
the basic 'naturalness' assumption imposed on the model, to avoid further fine tunings, 
is Mpi ~ M* ~ k. The geometry of Eq. (5) implies that the mass spectrum on the Planck 
brane, of the order of 10 15 TeV, can for kR c ~ 1 1 be exponentially 'warped' down to the 
TeV brane where SM particles live, by many orders of magnitude, namely, to the scale 
A^ = Mpiexp(— klzR c ) ~ 1 TeV. The appealing consequence is then that gravitational 
effects can occur in the reach of supercolliders, such as the LHC. Indeed, junction 
conditions at the brane y-positions imply, for the above value of kR c , the existence of 

(n) 

a tower of spin-2 graviton excitations, h^, with a specifically spaced mass spectrum 
M n = x n kexp (—knR c ) of order TeV (x n are here the roots of J\ (x n ) = 0). Their couplings 
to the SM particles are only 1/A^-suppressed (not 1 /Mpi): 



=^TeV = - 




1 Recently, an asymmetry Ace defined along these lines, in terms of differences between lepton and 
antilepton pseudorapidities, has been proposed for model identification in Ref. [11]. Also, the possibility 
of identifying the graviton excitation spin-2 from the azimuthal angular dependence of the graviton+jet 
inclusive production has been explored in [12]. 



and their signature may appear at LHC. In (6), T^ v is the energy-momentum tensor and 

h^l denotes the zero-mode, ordinary, graviton. 

The commonly chosen independent RS model parameters are the mass of the lowest 
graviton excitation, Mq = M\, and the 'universal', dimensionless, coupling c = fc/M P1 
(the scale A n is then a derived quantity). Theoretically 'natural' ranges of these param- 
eters are 0.01 < c < 0.1 and A n < 10 TeV [13]. With Y n = pM n x 2 n c 2 and p a number of 
the order of 0. 1, narrow graviton resonances are expected for such small values of c. 

Tevatron 95% CL limits on M G [14] range from M G > 560 GeV (c = 0.01) to 
Mq > 1.050 TeV (c = 0.1). Thus, due to the high allowed values of Mq, the discovery 
could be limited to the lightest M\ and the verification of the predicted mass pattern 
hardly feasible, so that the spin-2 determination through the angular analysis becomes a 
crucial test of the model. 

For dilepton production, the leading order subprocesses qq — >■ G —?■ and gg — > 
G — > give for the z-even distributions needed in Ace, with self-explaining notations 
[15,16]: 



d<7G = |(l+^ 2 )^ SM + ^-3^ 2 + 4z 4 )< + |(l-z 4 K 



dz 8 



8 



(7) 



and: 



^ce = e q 



SM iSM 



*CE 



2z* 5 + ^z*(l-z* 2 ) 



+ £ 



(8) 



In (8), £q, £g and £^ M are the fractions of G-events under the peak at Mr initiated by qq 
and gg processes, and the SM background, respectively. They are determined by overlaps 
of parton distribution functions and, obviously, £ q + Eg + e| M = 1. Strictly, Eqs. (7) 
and (8) are quite transparent in showing the characteristic z (and z*) dependencies for 
the spin-2 graviton, but hold in that form only for z cut = 1. This fact turns out to be 
numerically unimportant at the optimal value of z* where the Ace analysis is performed, 
and the final results take into account all phase space cuts foreseen by the experiment. 
Moreover, next-to-leading order terms in QCD [17] have for simplicity been included in 
the calculations through flat (in z) ^-factors, K ~ 1.3. 

For diphoton production (1), retaining only the leading order RS resonance exchange 
contributions to simplify the presentation, qq — >■ G — >■ yy and gg — > G — > yy, the ana- 
logues of Eqs. (7) and (8) with the same significance of the notations can be written as 
[18]: 

tiC7 ° 5 (l-z 4 )o? + J-(l + 6z 2 + z 4 )<7?, (9) 



dz 8 



32 



and 



A CE — fc o 



1 */c *4n 



H — z H — z H — z 
^8 4 8 



(10) 



It turns out that in this case the gg-initiated subprocess is the dominant, with a shape 
peaked at z = ±1 similar to the SM background but by far overwhelming it for the 
order TeV values of Mq we here are interested in. Clearly, the interest of the diphoton 



channel is that spin- 1 77 leaves only the spin-2 and spin-0 hypotheses, and that the 
ratio BR(G ->■ 77) /BR(G ->■ l+l~) ~ 2 is predicted [15]. 

A current field of interest is the extension of the simplest RS model described here 
to the case of SM fermions and gauge bosons also propagating in the full 5-D space, 
this could provide an approach also to other mass hierarchies such as, for example, the 
fermion mass hierarchy. A complicated phenomenology then emerges, in which both 
fermions and gauge bosons are accompanied by high (and very high) mass excitations, 
see for instance the review in Ref. [19]. The existence of the spin-2 graviton excitations 
remains, however, a general feature. 2 

Heavy neutral gauge bosons 

Turning to the spin-1 hypothesis, Z' models depend, besides the mass M = M Z ', on 
left- and right-handed couplings to quarks and leptons. In popular scenarios we refer 
to, generated by different extended electroweak gauge symmetries, those couplings are 
fixed theoretically, so that only M z r would remain as a free parameter. This is the case 
of the Z' x , Z'y, ZL, Z^ R , Z' ALR models, and the "sequential" Z' SSM model with the same 
couplings as the SM (details can be found in [2]). Current Tevatron lower limits (95% 
CL) on these Z' masses range from 878 GeV to 1.03 TeV, depending on the model [21]. 

Actually, the leading-order partonic subprocess qq — >■ Z' — >■ leads to the same 
form of the z-even angular distribution as the SM, therefore to the same Ace for all 
models: 

^ = |(i+^)[oT+of] (") 

and 

A* E =A™= 1 -z*(z* 2 + 3)-l. (12) 

Our AcE-based estimates for exclusion of the spin-1 Z' hypothesis will, accordingly, 
have a considerable degree of Z' model-independence. 

Sneutrino exchange 

The spin-0 character of the resonant subprocess d d — > v — > leading to a peak at 
M = My, manifests itself in a flat angular distribution [3]: 

d gV 3, , 2\_SM, 1 V /1 ox 

^=g(l+^K +2 G «> (13) 

A CE = ^ M A s ^ + <(2z*-l). (14) 

Besides M v , the cross section is proportional to the 7?-parity violating product X = 
{V) 2 Bi where 5/_is the sneutrino leptonic branching ratio and X' the relevant sneutrino 
coupling to the dd quarks. Current constraints on X are very loose (we may take the 
range 10~ 5 < X < 10 _1 ), and the 95% CL Tevatron lower limits on M v vary from 397 
GeV (X = 10~ 4 ) to 866 GeV (X = 10~ 2 ) [21]. 



2 The other one being represented by the spin-0, massive, radion needed to stabilize the compactification 
radius R c [20]. 




RESULTS FOR SPIN-2 IDENTIFICATION 

There are wide domains in Mr and coupling constant values allowed by current exper- 
imental limits where the scenarios described above predict the same peaks in M with 
the same numbers of events, so that they cannot be distinguished from each other on the 
basis of the event rates only. 

To proceed with the spin-identifying angular analysis, we suppose that a peak at 
M = Mr is discovered in process (1), and assume that it is consistent with a spin-2 RS 
graviton (in which case Mr = Mq). To evaluate the domain in the {Mq, c) plane, where 
the competitor spin-1 and spin-0 hypotheses with the same number of events under the 
M = Mq peak as the graviton hypothesis can be excluded, hence the spin-2 hypothesis 
can be established, we look at the 'distances' among models: 

A^CE = ^CE — ^CE an d A^CE = ^CE — ^CE- 0^) 

We adopt a simple-minded # 2 -like criterion, where the deviations (15) are compared 
to the expected statistical uncertainty dA^ E pertinent to the RS model (systematic 
uncertainties can also be included). The % 2 for the Z' and the sneutrino cases in (15) 
is defined as % 2 = \AA CE / 5Aq E \ 2 . The desired identification domain in (M<j,c) is 

determined by the condition % 2 > Xcl f° r ^ otn me ^ anc ^ me V models, with Xcl a 
critical number that specifies the confidence level (3.84 for 95% CL). Of course, Eq. (15) 
depends on the value of z*, but it turns out that the choice z* = 0.5 is 'optimal' in 
maximizing the sensitivity of the numerical analysis to the RS graviton. 

Figure 1 shows the 5-a discovery domains and the 95% identification domains 
for LHC luminosities of 10 and 100 fb 1 (Z = combined), and the 'theoretically 
favoured' restriction A n < 10 TeV taken into account. Specifically: the area at the left 
of the line "G" is the discovery domain for the lowest-lying RS graviton resonance; the 
area at the left of the "V" curve is the exclusion domain of the spin-1 hypothesis; finally, 
in the area at the left of the "5" line the spin-0 (and, as one can see, forcefully the spin-1) 
hypothesis can be excluded. Therefore, the identification domain where both spin-1 and 
spin-0 hypotheses can be excluded and the spin-2 established is the intersection of the 



TABLE 1. Discovery and Identification reach [TeV] 



Discovery Identification 





c = 0.01 


c = 0.1 


c = 0.01 


c = 0.1 


10 ftr 1 


1.7 TeV 


3.5 TeV 


1.1 TeV 


2.4 TeV 


100 ft- 1 


2.5 TeV 


4.6 TeV 


1.6 TeV 


3.2 TeV 




"5" area with the A n < 10 TeV one, i.e., above that curve. The allowed domain to the 
right of the dashed "oblique corrections" line is qualitatively determined by a fit to the 
oblique EW parameters [13, 22], but the condition dramatically restricting the discovery 
domains to the dashed areas, if literally applied, is the A n bound. This condition es- 
sentially forbids the creation of additional mass scale hierarchies in the model. Table 1 
summarizes the above results obtainable from the observation of RS graviton peaks in 
dilepton final states. 

Figure 2 shows some preliminary attempts to assess both discovery and identification 
reaches on the spin-2 RS graviton from the diphoton production channel in (1). Specifi- 
cally, the left panel shows, for 0.01 < c < 0. 1 the minimum number of events needed to 
identify by means of the asymmetry Ace a peak at diphoton invariant mass Mq as an RS 
resonance, against the spin-0 hypothesis (95% CL and LHC luminosity of 100 fb _1 ). 
Leading order QCD (therefore unit ^f-factors) has been used to obtain this figure. Also, 
the foreseen phase space experimental cuts for this channel have been taken into ac- 
count: p± > 40 GeV and |rj | < 2.4 for both photons, and reconstruction efficiency 80%. 
The right panel of Fig. 2 represents, in the same style as Fig. 1, the translation to the 
(Mq, c) RS parameter plane of the minimum number of diphoton events for RS graviton 
identification given in the left panel. This channel is still under study, in particular as 
regards the effect of the next-to-leading order QCD effects and the perspectives at lower 
LHC integrated luminosity. Nevertheless, the example in the figure indicates that dipho- 
ton events may have a comparable identification sensitivity to the spin-2 RS resonance, 
with the spin-1 hypothesis automatically excluded. 
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